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Energy  conservation  has  become  an  increasingly  important  issue  for  all  sectors,  particularly  in  industry. 
Therefore,  the  thermal  performance  of  insulation  systems  and  their  influence  on  heat  loss/gain  in  various 
applications  in  addition  to  economic  considerations  have  received  increased  attention  in  recent  years. 
In  this  study,  a  literature  review  of  papers  that  addressed  the  optimum  economic  thickness  of  the 
thermal  insulation  on  a  pipe  or  duct  with  different  geometries  used  in  various  industries  was  carried  out. 
The  studies  related  to  determining  the  critical  insulation  thickness  for  different  geometries  including 
circular  shapes  were  investigated.  The  heat  transfer  equations,  the  basic  results,  the  optimization 
procedures  and  the  economic  analysis  methods  used  in  the  studies  were  presented  comparatively. 
Additionally,  a  practical  application  example  based  on  optimizing  the  insulation  thickness  on  a  pipe  was 
performed,  and  the  effective  parameters  of  the  optimum  thickness  were  investigated. 

©  2013  Elsevier  Ltd.  All  rights  reserved. 
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1.  Introduction 

Thermal  insulation  systems  have  been  used  in  practice  for  many 
years  for  different  applications  and  purposes,  such  as  to  decrease  heat 
transfer  to/from  surfaces,  to  control  the  process  and  surface  tempera¬ 
tures,  to  avoid  the  condensation  problem,  and  to  provide  a  comfor¬ 
table  indoor  thermal  environment.  Increasing  concerns  regarding 
energy  efficiency,  climate  change  and  awareness  of  the  limited  energy 
resources,  the  use  of  a  proper  amount  of  thermal  insulation  for 
buildings  and  industrial  applications  has  gained  popularity.  When 
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examining  the  sectoral  distribution  of  the  energy  demand,  it  is  seen 
that  considerable  portions  of  the  global  energy  have  been  consumed 
in  the  residential  and  industrial  sectors,  accounting  for  approximately 
30%  and  40%,  respectively,  of  the  total  global  energy  [1-3].  Thermal 
insulation  is  primarily  used  to  limit  heat  loss/gain  from/to  surfaces 
under  operating  conditions  at  temperatures  above  or  below  ambient 
temperature,  i.e.,  to  provide  a  contribution  for  energy  conservation. 
Energy  conservation  is  a  major  concern  in  many  industrial  applica¬ 
tions.  The  primary  reasons  to  conserve  energy  include  maximizing 
the  return  on  investment  and  minimizing  the  life  cycle  cost  and  the 
emissions  associated  with  energy  consumption. 

In  addition  to  economic  considerations,  attention  has  focused 
on  improving  performance  and  thermal  efficiency  of  the  insulat¬ 
ion  systems  in  recent  years.  The  concept  of  an  economic  thermal 
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Nomenclature 

wp 

pump  power,  kj/s 

z 

capital  cost  rate,  US$/s 

Bi 

Biot  number  (/i  r0//<) 

<p 

exergy  efficiency  of  pipe  segment 

C 

cost,  US$/kg  or  US$/m3 

X 

insulation  thickness,  m 

C 

cost  rate,  US$/s 

rj 

heating  system  efficiency 

CRF 

capital  recovery  factor 

G 

Stefan-Boltzmann  constant,  W/m2I<4 

d 

discount  rate 

e 

surface  emissivity 

D 

diameter,  m 

6 

dimensionless  temperature  ( T-Tint)l(T[-T0 ) 

DP 

ratio  of  down  payment  to  initial  investment 

E 

exergy  rate,  kj/s 

Subscript 

h 

convective  heat  transfer  coefficient,  W/m2I< 

HDD 

heating  degree-days 

1 

pipe  inside  diameter 

Hu 

lower  heating  value  of  fuel,  kj/m3 

2 

pipe  outside  diameter 

i 

inflation  rate 

cr 

critic 

k 

thermal  conductivity,  W/mI< 

el 

electricity 

L 

pipe  length,  m 

en 

energy 

IT 

lifetime,  years 

/ 

fluid 

Ms 

ratio  of  first  year  miscellaneous  costs  to  initial  cost 

i 

inside 

N 

annual  operating  time,  hours 

ins 

insulation 

Pr 

Prandtl  Number 

int 

initial  condition 

PWF 

present  worth  factor 

max 

maximum 

Q 

heat  loss  rate  from  pipe,  kj/s 

min 

minimum 

r 

radius,  m 

0 

outside 

R 

thermal  resistance,  m2I</W 

opt 

optimum 

Rv 

ratio  of  the  resale  value  to  initial  investment 

p 

pump 

Ra 

Rayleigh  number 

s 

surface 

Re 

Reynolds  Number 

surr 

surrounding 

t 

wall  thickness,  m 

t 

total 

T 

temperature,  °C 

V 

velocity,  m/s;  volume,  m3 

insulation  thickness  considers  the  initial  cost  of  the  insulation  system 
plus  the  ongoing  value  of  the  energy  savings  over  the  expected  service 
lifetime  of  the  insulation.  Numerous  studies  on  thermal  insulation 
systems  for  different  applications,  including  buildings,  cold  stores, 
pipelines,  ducts,  tanks  and  other  equipment,  have  been  published  in 
the  literature.  The  majority  of  the  studies  that  are  related  to  optimizing 
the  thermal  insulation  thickness  focus  on  flat  surfaces,  such  as  building 
walls,  rather  than  cylindrical  geometries.  The  economic  insulation 
thickness  for  building  walls  depends  on  various  factors,  such  as 
building  type,  function,  wall  orientation,  construction  materials,  cli¬ 
matic  conditions,  insulation  properties  and  cost,  energy  type  and  cost 
and  efficiency  of  the  heating  or  air-conditioning  system  4-9].  In  these 
studies,  the  effect  of  mass  and  insulation  location  on  heating  and 
cooling  loads  were  analyzed  in  buildings  with  massive  exterior 
envelope  components  for  various  wall  configurations.  Considering 
the  heating  and/or  cooling  loads  of  buildings,  optimization  studies  on 
the  insulation  thickness  for  the  external  walls  of  buildings  have  been 
conducted  for  different  climate  zones  in  various  countries  including 
Qatar  [6],  China  [9],  Maldives  [10],  Turkey  [2,3,10-13],  Saudi  Arabia 
[14-16],  Malaysia  [17],  Palestine  [18]  and  Tunisia  49,20  .  A  suitable 
amount  of  thermal  insulation  in  the  building  envelope  can  result  in  a 
considerable  reduction  in  the  heating  and  cooling  energy  demands 
of  a  building  and  its  associated  C02  and  S02  emissions  into  the 
atmosphere.  In  addition  to  insulation,  several  investigations  have 
focused  on  optimizing  the  building  shape  to  minimize  the  energy 
demand  and  cost  [21-23].  All  of  these  studies  help  to  reduce  building 
energy  use  (annual  energy  requirements  for  heating  and  cooling)  and 
the  size  of  the  air-conditioning  and  heating  systems  in  buildings  and 
to  achieve  desirable  indoor  thermal  comfort  for  occupants. 

Compared  to  studies  on  flat  surfaces,  there  are  few  studies  related 
to  insulation  systems  on  pipes  and  ducts  despite  the  extensive 
applications  in  widely  diversified  fields.  Furthermore,  deriving  the 


heat  transfer  and  energy  cost  equations  and  taking  the  derivative  of 
the  objective  function  for  cylindrical  geometries  such  as  pipes  may  be 
considered  slightly  more  difficult.  In  most  studies,  optimum  insula¬ 
tion  thickness  computations  were  performed  based  mainly  on  the 
convective  and  radiative  heat  loss  from  a  pipe  or  duct  and  other 
parameters,  such  as  the  costs  of  the  insulation  material  and  energy, 
the  heating  system  efficiency,  the  lifetime  and  the  current  inflation 
and  discount  rates.  The  heat  loss  from  a  pipe  (or  the  energy 
requirement  to  heat  the  fluid  in  a  pipe)  is  the  main  input  required 
to  analyze  the  optimum  insulation  thickness.  Heat  loss  occurs  by 
conduction,  convection  and  thermal  radiation.  Most  of  the  studies 
considered  only  convective  [24-31  or  radiative  heat  transfer  modes 
[32]  while  several  studies  considered  both  modes  [33-36].  Moreover, 
Kecebas  et  al.  [37]  and  Basogul  and  Kecebas  [38]  used  the  degree-days 
(DD)  for  estimating  the  heating  energy  requirement. 

In  this  paper,  studies  related  to  determining  the  optimum  thermal 
insulation  thickness  for  pipes  and  ducts  were  reviewed  first.  The 
optimization  procedure  was  introduced,  and  the  operating  conditions 
and  parameters  used  in  these  studies  were  listed.  Then,  studies 
determining  the  optimum  insulation  thickness  on  pipes  and  ducts 
with  various  geometries  on  reducing  convective  and  radiative  heat 
transfer  were  investigated  and  presented.  After  examining  the  eco¬ 
nomic  analysis  methods  used  to  obtain  an  economic  insulation 
thickness,  a  simple  and  practical  application  for  determining  the 
optimum  insulation  thickness  for  pipes  was  conducted. 

2.  Determining  the  economic  insulation  thickness 

2.1.  Optimization  procedure 

To  minimize  the  energy  and  insulation  costs  in  addition  to 
reducing  the  heat  loss  to  the  surroundings,  the  thickness  of  the 
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insulation  material  needs  to  be  optimized.  The  motivation  for 
including  pipe  insulation  is  often  to  minimize  the  total  cost,  which 
includes  the  cost  of  the  insulation,  its  installation  and  maintenance 
as  well  as  the  cost  of  the  energy  loss  by  convective  or  radiative 
heat  transfer.  Heat  transfer  principles  and  cost  information  are 
used  in  defining  the  overall  cost  function  to  be  minimized  [33]. 
The  economic  insulation  thickness  for  a  pipe  is  a  function  of  a 
large  number  of  parameters,  such  as  pipe  size,  cost,  conductivities 
of  the  pipe  and  the  insulation  material,  operating  and  ambient 
temperatures,  heat  transfer  coefficients  at  the  inside  and  outside 
of  the  pipe,  economic  parameters  and  annual  operation  hours 
[25,29,38,39].  The  optimum  economic  thickness  is  generally 
accepted  as  the  value  that  provides  the  minimum  total  life  cycle 
cost,  as  illustrated  in  Fig.  1.  As  the  thermal  insulation  thickness  on 
a  pipe  increases,  the  heating  and  cooling  transmission  loads  and 
their  costs  (i.e.,  energy  costs)  decrease.  The  transmission  loads  are 
used  as  the  input  data  for  an  economic  model  to  determine  the 
variation  in  the  cost  of  the  insulation  plus  the  present  value  of  the 
energy  consumption,  which  is  considered  lost  energy,  over  the 
lifetime  of  the  system.  Conversely,  as  shown  in  Fig.  1,  the  insula¬ 
tion  cost  increases  with  the  amount  of  material  used.  In  contrast  to 
flat  surfaces,  the  variation  in  the  amount  of  insulation  used  on  a 
pipe  versus  the  insulation  thickness  is  not  linear,  but  rather,  the 
increase  is  parabolic. 


2.2.  Optimum  thermal  insulation  thickness  for  pipes  and  ducts 

Economic  evaluations  regarding  thermal  insulation  in  district 
heating  pipelines  were  investigated  by  Kecebas  et  al.  [37]  using  the 
heating  degree-days  (HDD)  indicator  in  Afyonkarahisar/Turkey. 
In  the  calculations,  coal,  natural  gas,  fuel-oil  and  geothermal  energy 
were  used  as  the  energy  sources,  and  rock  wool  was  used  as  the 
insulation  material  for  nominal  pipe  sizes  (ranging  from  50  mm  to 
250  mm).  Eqs.  (1)  and  (2)  were  used  for  the  convective  heat  transfer 
coefficients  inside  and  outside  the  insulated  pipe,  respectively.  The 
results  showed  that  the  highest  values  of  the  optimum  insulation 
thickness  and  the  energy  savings  were  reached  for  a  nominal  pipe 
size  of  250  mm  with  an  insulation  thickness  of  0.228  m.  Similarly,  the 
environmental  evaluations  of  the  thermal  insulation  in  the  district 
heating  pipeline  were  investigated  by  Basogul  and  Kecebas  [38]. 
Using  thermal  insulation  on  pipelines  helps  to  reduce  the  heating  or 
cooling  energy  requirements  of  transmitted  fluid  and  its  associated 
greenhouse  gas  emissions  into  the  atmosphere.  The  reductions  in  the 
emissions  of  C02,  CO  and  S02  with  different  insulation  thickness  for 
nominal  pipe  sizes  ranging  from  50  mm  to  200  mm  were  examined 
under  the  same  conditions  as  those  used  in  Kecebas  et  al.  [37].  It  was 
found  that  when  the  optimum  thermal  insulation  was  applied  to  a 


heating  pipeline,  the  emission  of  C02  was  decreased  by  21%. 

/ii  =  0.023Re08Pr°4(fe//D1)  (1) 
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where  Dt  and  Dins  are  the  inside  pipe  and  insulation  diameters, 
kf  is  the  thermal  conductivity  of  the  fluid,  TSi0  is  the  mean  outer  surface 
temperature  of  the  insulation,  T0  is  the  outside  air  temperature  and  V 
is  the  outside  air  velocity,  which  was  assumed  to  be  0.2  m/s. 

Ozdemir  and  Parmaksizoglu  [24]  focused  on  determining  the 
optimum  insulation  thickness  for  pipes  used  in  low  temperature 
heating  systems  in  buildings.  The  variation  in  the  annual  net 
savings  with  the  insulation  thickness  was  investigated  using  the 
life  cycle  cost  analysis,  but  the  effects  of  the  pipe  diameter,  the 
operation  period,  the  operation  temperature  of  water  and  the 
properties  of  the  insulation  material  were  not  considered  in  that 
study.  The  convective  heat  transfer  coefficient  inside  the  pipe  was 
assumed  to  be  constant  and  equal  to  4000  W/m2K.  As  opposed  to 
the  study  by  Kecebas  et  al.  [37],  the  following  equations  were  used 
for  the  convective  heat  transfer  coefficient  outside  the  insulated 
pipe  for  laminar  and  turbulent  flow,  respectively. 
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where  AT  is  the  temperature  difference  between  the  outer  surface 
of  insulation  ( TinStS )  and  the  outside  air  (T0)  in  °C,  rins  is  the 
insulated  pipe  radius  in  m  and  V  is  the  outside  air  velocity  in  m/s. 

In  the  study  by  Li  and  Chow  [25],  the  thermal  insulation  and 
heating  cable  systems  were  used  together  to  protect  the  water 
pipes  from  freezing  in  cold  regions  such  as  Beijing/China.  In  the 
calculations,  the  variation  in  the  optimum  insulation  thickness 
was  determined  depending  on  the  average  outside  air  tempera¬ 
ture,  the  outer  pipe  diameter,  the  lifetime,  the  thermal  conductiv¬ 
ity  of  the  insulation  material  and  the  insulation  cost.  An  empirical 
formula  given  by  Eq.  (5)  was  used  to  estimate  the  convective 
heat  transfer  coefficient  as  a  function  of  the  wind  speed  of  the 
environment.  The  average  outdoor  air  temperature  and  the  mean 
wind  speed  were  assumed  to  be  - 10  °C  and  4.8  m/s,  respectively. 
It  was  found  that  increasing  the  lifetime  of  the  system  and 
decreasing  the  outside  air  temperature  result  in  a  thicker  optimum 
insulation  thickness.  Furthermore,  the  optimum  insulation  thick¬ 
ness  was  inversely  proportional  to  the  thermal  conductivity  cost  of 
the  insulation. 


h0  =  11.73  +  l/1/2  (5) 

For  cold  applications,  the  C02  emissions  and  the  energy  savings 
potential  obtained  by  using  correct  pipe  insulation  in  space 
heating  and  domestic  hot  water  pipelines  were  investigated  by 
Chmielarski  [40  .  In  the  study,  the  ambient  air  temperature,  the 
chilled  water  supply  temperature  and  the  air  conditioning  period 
were  assumed  to  be  26  °C,  7  °C  and  6  months,  respectively. 
Assuming  typical  European  basic  conditions,  the  optimal  insula¬ 
tion  thickness  varied  from  20  mm  to  30  mm.  Information  regard¬ 
ing  the  studies  is  summarized  in  Table  1. 

The  outside  heat  transfer  coefficient  equations  for  horizontal 
and  vertical  steel  pipes  as  a  function  of  the  surface  temperature  of 
the  steel  pipe  (Ts)  and  the  ambient  air  temperature  (T0)  are  given, 
respectively,  as  follows  [30]: 

h0  =  3.2  +  0.05(TS-  T0) 


(6) 


Table  1 

Summary  information  about  the  studies  related  to  the  optimization. 


Study 

Parameters 

Economic  parameters 

Operating  conditions 

Notes 

Kecebas  et  al.  [37] 

D2  =  50-250  mm,  kpipe=  16.2  W/mK, 
kins  =  0.040  W/mK  (rock  wool), 

Cins= 95  US$/m3 

P!-P2  method  was  used.  LT=  10 
years, 

interest  and  inflation  rates  are  4% 
and  5%  respectively. 

T0=15  °C,  HDD= 2828,  energy  sources  are:  coal,  natural  gas, 
fuel-oil  and  geothermal  energy. 

Basogul  and  Kecebas  [38] 

D2  =  50-200  mm,  kpipe= 54  W/mK, 

P1-P2  method  was  used.  LT=  10 

HDD= 2828,  energy  sources  are:  coal, 

Environmental  impacts  were  also 

kfns= 0.040  W/mK  (rock  wool), 

Cins— 95  US$/m3 

years, 

interest  and  inflation  rates  are  4% 
and  5%  respectively. 

natural  gas,  fuel-oil,  LPG  and  geothermal. 

investigated. 

Ozd  emir  and 

D2  =  108  mm,  kpipe=40  W/mK, 

LT=  15  years,  interest  rate  is  5%. 

Tf=  80  °C,  T0  =  20  °C,  h,  =  4000  W/m2K, 

Fiberglass  and  natural  gas  were  assumed  as 

Parmaksizoglu  [24] 

kfns=0.04  W/mK  (fiberglass) 

annual  operation  time  is  5040  h. 

insulation  and  energy  source.  The  heating 
system  efficiency  was  assumed  to  be  0.9. 

Li  and  Chow  [25] 

D2  =  20-200  mm, 

LT=  5-10  years,  interest  and 

T0=  — 10  °C,  Ts„o=  -26  °C,  h0  =  13.9  W/m2I< 

The  total  cost  includes  the  heating  cable 

kins =0.024  W/mK 

inflation 

rates  are  4%  and  5%  respectively. 

(for  V=4.8  m/s),  3000  h. 

cost 

together  with  the  operating  energy  and 
insulation  costs. 

Kecebas  [42] 

D2  =  50  mm,  kpipe = 54  W / ml( 

(stainless  steel), 

kins= 0.040  W/mK  (rock  wool),  C,-ns= 95  US 
$/m3 

LT=  10  years,  interest  and  inflation 
rates  are  4%  and  5%  respectively. 

Tf  =90  °C 

Energy  sources  are  natural  gas,  coal 
and  fuel-oil.  The  optimization  method  is 
based  on  exergo-economic  analysis. 

Chmielarski  [40] 

D2=40-115  mm 

Energy  cost  increase  is  5%,  interest 
rate  is  4%. 

T0= 26  °C,  Tf=  7  °C,  operating  period 

6  months,  energy  efficiency  ratio  2.5. 

Ozturk  et  al.  [36] 

100  <  D2  <  600  mm,  0  <  xins  <160  mm, 
kPipe  =  54  W/mK,  kins= 0.045  W/mK 

Interest  rate  is  6%,  LT  =20  years. 

Fluid  is  water,  Tf=  120  °C,  To  =  10  °C,  energy 

source  is  natural  gas,  annual  operation  time  is  3200  h. 

Soponpongpipat  et  al.  [26] 

D2  =  0.5  m,  kdUCt—  60.5  W/mK, 

P1-P2  method  was  used,  interest 

T/-T0  =  13  °C,  fluid  is  air,  average  COP  of  air  conditioning  system 

Rubber  and  fiberglass  were  assumed  as 

kins  =  0.035  W/mK  (rubber), 
kim— 0.045  W/mK  (fiberglass) 

rate  is  14.25%,  inflation  is  2.1%, 

LT=  10  years. 

is  2.81,  total  operation  time  23040  h. 

main  and  auxiliary  insulation  materials, 
respectively. 

Zaki  and  Al-Turki  [29] 

D2  =  10-27.3  cm,  insulation  materials: 
rock  wool  and  calcium  silicate. 

Interest  rate  was  not  given. 

8760  h,  ho  =  10  W/m2K,  Tf=  170-400  °C. 

Superheated  systems  were  studied,  such 
as  steam,  furfural,  crude  oil  and  300- 
distillate  oil. 

Soylemez  [43] 

Cduct=10  $/kg,  Cf=  0.075  $/kWh  (electricity) 

P1-P2  method  was  used,  energy 
price  and  interest  rate  are  0.1, 

LT=  20  years. 

Annual  operating  hours  is  1600  h. 

A  present  sample  problem  was  used  for 
comparison. 

Kalyon  and  Sahin  [27] 

kinS =0.036  W/mK,  kpipe=30  W/mK, 
r PiPe =0.0127  m,  tpipe= 0.002  m 

Economic  analysis  was  not  used. 

ho  =  10  W/m2I<,  To  =  0  °C,  Tf=  77  °C. 

Only  convective  heat  transfer  was 
considered. 

Sahin  [32] 

ki„s =0.036  W/mK,  kpipe=30  W/mK, 
r pipe  =  0.0127  m,  tpipe= 0.002  m 

Economic  analysis  was  not  used. 

Ta=  -273  °C,  Tf=  77  °C,  ht= 122  W/m2K,  s  =  0.80. 

Only  radiative  heat  transfer  was 
considered. 

Sahin  and  Kalyon  [34] 

kfns  =  0.01  ■ -0.05  W/mK, 

kp^e  =  30  W/mK,  rp!pe= 0.0127  m 

Economic  analysis  was  not  used. 

ho  =  10  W/m2K,  hi=  102.6  W/m2K,  To  =  0  °C,  Tf=  77  °C. 

Convective  and  radiative  heat  transfer  were 
considered. 

Bahadori  and  Vuthaluru  [44] 

D2  =  0.02-0.5  m, 
k,-ns= 0.02-0.09  W/mK 

Economic  analysis  was  not  used. 

Ts=  100-700  °C,  Ta= 20  °C 

The  average  absolute  deviation  is 
approximately  2.12%  depending  on  surface 
temperature  and  diameter. 

Bahadori  and  Vuthaluru  [45] 

100  <  D2  <  2400  mm 

Economic  analysis  was  not  used. 

Ts-T0  <  250  °C 

The  average  absolute  deviation  is 
approximately  3.25%. 

Sahin  et  al.  [46] 

D2  =  0.1  m,  kPipe =0.144  W/mK 

The  costs  of  irreversibilities  due  to 
the  heat  transfer  and  pressure 
drop  were  computed. 

Ts=350  K,  Ta= 300  K,  L=100  m,  kf=  2.25  W/mK 
(engine  oil) 
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fto=3.4  +  0.09(rs- T0)  (7) 

Ozturk  et  al.  [36]  compared  four  optimization  methods  for  the 
design  of  hot  water  piping  systems.  The  aim  of  the  first  method  is 
to  optimize  both  the  pipe  diameter  and  the  insulation  thickness. 
First,  the  pipe  diameter  is  optimized  on  the  basis  of  the  pipe 
investment  cost  and  the  operation  cost  due  to  the  friction  of 
the  fluid  in  the  pipe.  Then,  the  insulation  thickness  on  a  pipe  is 
optimized  by  considering  energy  and  insulation  costs.  The  opera¬ 
tion  cost  is  related  to  the  power  required  to  pump  the  fluid  along 
the  piping.  The  insulation  cost  is  a  function  of  the  thickness  of  the 
insulation  that  is  applied  on  the  pipe.  The  energy  cost  is  a  concept 
that  is  related  to  the  energy  source  type  and  the  heat  loss  per 
length  of  pipe.  In  Ozturk  et  al.  [36],  natural  gas  was  used  as  an 
energy  source.  The  objective  functions  to  optimize  the  pipe 
diameter  and  the  insulation  thickness  are  given,  respectively,  as 
follows: 


pipe  —  Zpipe  “F  (^e/Wp) 

(8) 

ins  =  Zins+iChQ) 

0) 

where  C  is  the  cost  rate  in  US$/s,  Z  is  the  capital  cost  rate  in  US$/s, 
Cei  and  Ch  are  the  specific  cost  for  electricity  and  the  unit  energy 
cost  of  hot  water  in  US$/kJ,  respectively,  Wpis  the  pump  power  in 
kj/s  and  Q.  is  the  rate  of  heat  loss  from  the  pipe  in  kj/s.  The  capital 
cost  rates  for  the  pipe  and  the  insulation  can  be  expressed  by 
Eqs.  (10)  and  (11)  as  follows: 

Zpipe  —  (FRF  4-  <;)Cpipe  (10) 

Zins=(CRF  +  g)Cins  (11) 

where  CRF  is  the  capital  recovery  factor,  which  is  dependent  on  the 
annual  interest  rate  and  the  lifetime  (LI),  q  is  a  coefficient  that 
accounts  for  part  of  the  fixed  operation  and  the  maintenance  cost, 
which  was  assumed  to  be  0.01,  and  Cpipe  and  Cins  are  the  pipe  and 
insulation  costs,  respectively. 

The  second  method  considers  not  only  the  insulation  but  also 
the  pipe  system,  and  it  is  based  on  the  minimization  of  the  total 
cost.  The  total  cost  includes  the  pipe  and  insulation  investment 
costs  and  the  operating  costs  due  to  pressure  and  heat  losses, 
which  is  given  by 

Q  =  (Zpipe  +Zins)(CeiWp  +  C/7Qj  (12) 

The  third  and  fourth  methods  mentioned  in  Ozturk  et  al.  [36] 
involve  an  exergy  analysis  for  the  water  stream  in  a  pipe.  The  third 
method  is  focused  on  the  minimization  of  entropy  generation  or 
the  maximization  of  exergy  efficiency  without  cost  effects.  The 
objective  of  the  method  is  to  increase  the  coefficient  of  thermal 
performance  (CTP),  and  which  is  given  as  follows: 

CTP=rr^  (13) 

0  -<p) 

where  cp  is  defined  as  the  exergy  efficiency  of  the  pipe  segment, 
which  is  the  ratio  of  the  exergy  flow  at  the  pipe  outlet  to  the 
exergy  flow  at  the  pipe  inlet  plus  the  work  supplied  for  hot  water 
transportation.  The  pipe  diameter  and  the  insulation  thickness  are 
determined  considering  the  maximum  CTP. 

The  fourth  and  most  recent  method  mentioned  in  Ozturk  et  al. 
[36]  considers  the  costs  of  exergy  destruction  due  to  the  friction 
loss  and  the  exergy  loss  from  heat  loss  in  addition  to  the  pipe  and 
insulation  investment  costs.  The  objective  function  used  to  deter¬ 
mine  the  economic  sizes  of  the  pipe  and  the  insulation  is  given  by 
the  following: 

Q  =  ( Zpipe  +Zim)(C  eldest  +  FeEi0SS )  (14) 

where  Cel  and  Ce  are  the  specific  exergy  costs  for  electricity  and  hot 
water,  respectively,  Edest  and  £/oss  are  the  exergy  destruction  rate 


due  to  friction  in  the  pipe  and  the  exergy  loss  rate  due  to  heat  loss 
from  the  pipe  to  the  environment,  respectively.  The  details  of  the 
equations  can  be  found  in  Ozturk  et  al.  [36]  and  Karabay  [41]. 
Using  these  equations,  Ozturk  et  al.  [36]  compared  the  results  of 
these  methods  when  using  the  same  criteria  in  a  case  of  a  pipe 
system.  Several  parameters  used  in  that  study  and  their  numerical 
values  are  given  in  Table  1.  As  a  conclusion,  in  the  study,  the  fourth 
method  was  recommended  to  be  used  in  the  design  studies 
because  the  exergy  and  cost  parameters  were  considered.  Simi¬ 
larly,  in  the  study  of  Kecebas  [42],  the  exergoeconomic  approach 
was  used  for  optimizing  the  pipe  insulation  thickness.  The  exer- 
getic  costs  of  insulation  and  fuels  (coal,  natural  gas  and  fuel-oil) 
were  calculated  based  on  life-cycle  cost  analysis.  The  parameters 
considering  in  the  analysis  were  excess  air  coefficient,  the  inlet 
temperature  of  fuel,  the  temperature  of  combustion  chamber  and 
the  temperature  of  stack  gases.  It  was  stated  in  the  study  that  the 
optimum  insulation  thickness  for  the  exergoeconomic  optimiza¬ 
tion  was  higher  than  that  of  energoeconomic  optimization. 

Soponpongpipat  et  al.  [26]  studied  the  economic  thickness  of  a 
double-layer  insulation  composed  of  the  rubber  (/<= 0.035  W/mI<) 
and  fiber  glass  (/<= 0.045  W/mK)  for  an  air  conditioning  duct. 
The  effects  of  the  convective  heat  transfer  coefficients  of  the  inside 
and  outside  of  a  circular  galvanized  steel  duct  on  the  optimum 
thickness  of  these  insulators  were  investigated  when  the  heat 
transfer  coefficients  varied  from  6  W/m2I<  to  22  W/m2I<.  It  was 
found  that  the  variation  in  the  inside  and  outside  heat  transfer 
coefficients  did  not  affect  the  optimum  thickness,  but  the  net 
savings  increased  with  increasing  the  heat  transfer  coefficients. 

Zaki  and  Al-Turki  [29]  presented  a  model  for  calculating  the 
optimum  thickness  of  an  insulation  system  in  which  the  compo¬ 
site  materials  have  different  properties  and  costs  for  pipes  (dia¬ 
meter  of  10  cm  -27.3  cm)  with  a  flow  of  superheated  steam. 
Mineral  rock  wool  and  calcium  silicate  were  assumed  to  be  the 
insulation  materials. 

In  addition  to  optimization  of  insulation  thickness,  Soylemez 
[43]  studied  the  optimum  duct  sizing  for  heating,  ventilation  and 
air  conditioning  (HVAC)  systems  with  round  and  rectangular 
shapes.  It  was  stated  that  the  airflow  velocity  in  an  HVAC  channel 
is  directly  related  to  the  operational  cost  of  the  HVAC  system  along 
with  the  initial  system  cost.  The  derivative  of  the  total  life  cycle 
cost,  which  is  the  sum  of  the  initial  cost  and  the  operating  (energy) 
cost,  was  taken  with  respect  to  the  pressure  drop  per  unit  channel 
length.  The  optimum  duct  diameter  was  determined  to  obtain  the 
minimum  pressure  drop.  In  the  optimization,  the  P1-P2  economic 
method  was  used  and  the  technical  lifetime  of  the  duct  was 
assumed  to  be  20  years.  Sahin  et  al.  [46]  numerically  investigated 
the  effect  of  fouling  thickness,  which  leads  to  decreased  heat 
transfer  and  increased  frictional  pressure  drop,  on  the  operational 
cost  related  to  irreversibilities  in  pipe  flow.  The  numerical  model 
developed  was  based  on  entropy  generation  due  to  heat  transfer 
and  viscous  friction. 

All  of  the  studies  mentioned  above  are  given  in  Table  1,  which 
describes  the  parameters,  assumptions  and  operating  conditions 
of  the  systems. 

2.3.  Studies  not  including  an  economic  analysis 

Heat  transfer  modes  depend  mainly  on  the  surface  tempera¬ 
tures  and  the  geometry  of  the  system  being  considered.  Pipes  or 
ducts  surrounded  by  gases  may  transfer  heat  via  conduction, 
convection,  and/or  thermal  radiation.  For  air,  conduction  is  usually 
dominated  by  other  modes  and  thus,  may  be  neglected  [47]. 
Kalyon  and  Sahin  [27]  studied  the  optimum  insulation  thickness 
of  a  circular  pipe  subjected  to  convective  heat  transfer.  The 
objective  is  to  determine  the  insulation  thickness  for  a  limited 
amount  of  insulation  material,  which  minimizes  the  heat  loss. 


0.  Kaynakli  /  Renewable  and  Sustainable  Energy  Reviews  30  (2014)  184-194 


189 


For  this  reason,  the  mentioned  study  did  not  use  any  economic 
methods.  They  found  that  a  uniform  insulation  layer  on  the  surface 
of  a  pipe  provides  optimum  results  for  the  case  of  a  convective 
heat  transfer  mode.  Sahin  [32  numerically  optimized  the  variation 
in  the  thermal  insulation  thickness  of  a  pipe  for  space  applications 
to  minimize  the  radiative  heat  loss  to  the  ambient.  The  thickness 
of  the  insulation  was  assumed  to  be  linearly  varying  along  the  pipe 
because  this  scenario  is  easy  to  implement  in  practice.  The  study 
aimed  to  determine  the  suitable  slope  of  the  insulation  thickness 
along  the  pipe,  which  maximizes  the  fluid  outlet  temperature 
(thus  minimizes  the  heat  loss).  It  was  found  that  a  linearly 
decreasing  insulation  thickness  with  a  minimum  slope  along  the 
pipe  provides  the  best  insulation  under  the  radiation  heat  transfer 
condition. 

As  reported  in  Sahin  and  Kalyon  34],  the  weight  and  volume  of 
the  insulation  material  are  required  to  be  minimal  in  some  cases, 
such  as  in  aerospace  applications.  For  this  reason,  they  studied  the 
insulation  thickness  variation  over  a  pipe  transporting  a  high 
temperature  fluid  to  maintain  a  uniform  outer  surface  tempera¬ 
ture.  The  heat  transfer  on  the  outer  surface  of  the  insulation  was 
considered  to  be  a  combination  of  convection  and  radiation.  In  the 
study,  the  variation  in  the  insulation  thickness  over  the  pipe  with  a 
dimensionless  axial  distance  was  calculated.  The  analytical  solu¬ 
tion  showed  that  the  insulation  thickness  variation  providing  the 
uniform  surface  temperature  is  independent  of  the  surrounding 
thermal  conditions  (i.e.,  the  convection  and  radiation  heat  transfer 
coefficients)  provided  that  these  conditions  are  uniform  over  the 
surface.  Moreover,  the  insulation  thickness  variation  along  a  pipe 
was  determined  in  an  exponential  form  (decreasing  in  an  expo¬ 
nential  function);  however,  it  was  very  close  to  a  linear  variation. 

Bahadori  and  Vuthaluru  [44]  formulated  a  simple  correlation 
to  estimate  the  economic  thickness  of  the  thermal  insulation  for 
pipes.  The  proposed  correlation  depends  on  the  steel  pipe  dia¬ 
meter  and  the  thermal  conductivity  of  the  insulation  and  covers  a 
pipeline  diameter  and  a  surface  temperature  of  up  to  0.5  m  and 
700  °C,  respectively.  Similarly,  Bahadori  and  Vuthaluru  [45] 
developed  a  simple  method  to  estimate  the  thermal  insulation 
thickness  required  to  arrive  at  a  desired  heat  flow  or  surface 
temperature  for  flat  surfaces,  ducts  and  pipes.  The  proposed 
method  considers  temperature  differences  between  the  ambient 
and  outside  temperatures  up  to  250  °C,  the  temperature  drop 
through  the  insulation  up  to  1000  °C  and  the  outside  diameter  of  a 
pipe  up  to  2400  mm.  It  was  noted  in  the  study  that  the  results  of 
the  model  are  in  good  agreement  with  the  reported  data  for  a 
wide  range  of  conditions,  and  the  average  absolute  deviation  was 
found  to  be  3.25%  for  a  temperature  drop  between  25  °C  and 
1000  °C  through  the  insulation. 

Alawadhi  31]  investigated  the  use  of  a  phase  change  material 
(PCM)  as  thermal  insulation  for  a  pipe.  The  effectiveness  of  the 


PCM  insulation  was  evaluated  by  comparing  its  thermal  perfor¬ 
mance  with  that  of  insulation  without  a  phase  change.  N- 
octadecane  (Ci8H38)  paraffin  was  selected  as  the  PCM.  Both  time- 
dependent  and  time-independent  boundary  conditions  were 
examined.  At  the  outer  surface,  the  radiation  effect  was  not 
considered  and  the  free  convection  boundary  condition  (0O  =  0, 
Bi  =  0.5)  was  assumed.  The  results  indicated  that  for  the  time- 
independent  case,  the  PCM  insulation  is  effective  in  reducing  the 
heat  loss  from  the  pipe  if  the  Rayleigh  number  is  low 
(Ra  =  5  x  106);  but  for  Ra  =  l  x  107,  the  PCM  insulation  shows  a 
weak  thermal  performance.  For  the  time-dependent  case,  the  heat 
loss  is  effectively  reduced  with  the  PCM  insulation  for  the 
examined  range  of  oscillating  periods. 


2.4.  Energy,  insulation  and  total  cost  equations 

The  total  thermal  resistance  for  the  circular  pipe  illustrated  in 
Fig.  2  can  be  written  as  follows: 

zE_(  '  +  (15) 

^7rKpipgL  ^jrVijjgLrio  J 

where  and  r2  are  the  inner  and  outer  radii  of  the  pipe, 
respectively,  hi  and  hQ  are  the  convective  heat  transfer  coefficients 
of  the  inner  and  outer  surfaces  of  the  pipe,  respectively,  L  is  the 
length  of  the  pipe,  rins  is  the  radius  of  the  insulated  pipe,  kpipe  and 
kins  are  the  thermal  conductivities  of  the  pipe  and  insulation 
materials,  respectively.  The  rate  of  the  heat  loss  from  the  pipe 
per  unit  length  can  be  calculated  as  follows  [27,32,33,35,37,48]: 

Q.  = _ 2 nAT _ =  (Tf-T0) 

L  (1  /v\  h{)  +  (/n(r2/r-i ) /kpipe) 3-  iio(r  [ns  /  v2)  /  k\ns)  +(l/rl-rjS/i0)  XR 

(16) 

where  7}  and  T0  are  the  temperatures  of  the  fluid  in  the  pipe  and 
the  outside  air. 

The  energy  cost  due  to  the  heat  loss  rate  is  given  as  follows: 

Cen2x3600ATN 

en,t  ~  ((1  /ri  hi)  +  ( In(r2 /rt )/kpipe)  +  (In(rins/r2)/kins)  +  (1  /rinsh0))Hu  n 

(17) 

where  Cen  is  the  unit  cost  of  the  energy  source  in  US$/kg,  US$/m3 
or  US$/kWh,  depending  on  the  energy  source  type,  N  is  the  annual 
operating  time  in  hours,  Hu  is  the  heating  value  of  the  fuel  in  kj/kg 
or  kj/m3,  depending  on  the  energy  source  type  and  rj  is  the  heating 
system  efficiency,  which  was  generally  assumed  to  be  0.90  or  0.93 
for  natural  gas  usage  [2,11,24,36-38,49,50],  0.65,  0.70  or  0.77  for 
coal  [2,11-13,37,38,51],  0.80  for  fuel  oil  [2,11,30,37,38,50]  and  0.99 
for  electricity  2,11  . 


ins 


Fig.  2.  An  insulated  pipe. 
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The  amount  of  used  insulation  per  unit  length  of  pipe  and  its 
cost  are,  respectively,  given  as  follows: 


ins 


=  xOins  ~  W 


Cins,t  —  Q  ns(Yins/L)  —  Qns^C^/ns  C?) 


(18) 

(19) 


where  Vins  is  the  volume  of  insulation  material  used  in  m* * 3,  Cins  is 
the  unit  cost  of  the  insulation  in  US$/m3  and  Cins  t  is  the  cost  of 
insulation  per  unit  length  of  pipe  in  US$/m. 

After  calculating  the  energy  and  insulation  costs,  the  present 
worth  of  the  total  life  cycle  cost  associated  with  the  pipe  insulation 
thickness  can  be  expressed  as  follows: 


Cen2^3600  AT  N  PWF 

((1  /r\  hi)  +  (In(r2  Ai )/kpipe)  +  (Jn(rins/r2)/kfns)  +  (1  /rinsh0))Hu  n 


+  C ins  Mr  ins 


(20) 


where  PWF  is  the  present  worth  factor,  the  formulation  of  which  is 
given  in  the  Economic  Analysis  section.  The  optimum  insulation 
thickness  is  obtained  by  minimizing  the  total  cost.  The  minimum 
point  can  be  determined  by  setting  the  first  derivative  of  the  total 
cost  function  to  zero  with  respect  to  the  radius  of  the  insulated 
pipe.  Consequently,  the  partial  derivative  of  Ct  with  respect  to  rins 
can  be  computed  ( SCt/5rins ).  Equating  the  first  derivative  function 
to  zero  yields  the  value  of  rinSt0pt,  which  is  determined  as  follows: 

o  =  (h0rins,opt  -  k/ns)1/2/(nns,opt)3/2((3600Cen(/irN)PWF)/(Hi////<irishoCins))1/2 
—  ((1  /f\  hi)- f-  (In(r2/r^  )/kpipe))  —  (Oti(rins  0Pf  /  r2)  /  k(ns)  3-  (1  /^ms,optho)) 

(21) 

Hence,  the  optimum  insulation  thickness  for  a  pipe  can  be 
written  as  follows: 

Xopt  —  fins, opt  —  l~2  (22) 


3.  Studies  on  heat  transfer  characteristics  in  ducts  and  pipes 

Several  investigations  have  focused  on  the  heat  transfer  char¬ 
acteristics  in  various  geometries  (circular,  triangular,  rectangular, 
oval,  polygonal  ducts  and  double  pipes)  rather  than  optimizing  the 
insulation  thickness  by  considering  the  economic  parameters 
[28,33,35,48,52-55  .  Chen  and  Yang  [54]  investigated  the  heat 
transfer  rate  on  the  insulation  layer  of  a  double  circular  pipe  and 
the  temperature  distribution  in  the  double  pipe.  Wong  et  al.  [55] 
analyzed  the  two  dimensional  steady  state  heat  transfer  charac¬ 
teristics  of  an  insulated  triangular  duct  by  using  the  two  models, 
which  are  the  plane  wedge  thermal  resistance  and  plate  thermal 
resistance  models.  The  critical  and  neutral  thicknesses  as  well  as 
heat  transfer  characteristics  for  an  insulated  pipe  have  been 
interesting  subjects  to  researchers.  The  critical  and  neutral  heat 
transfer  characteristics  were  investigated  by  Chou  [52]  and  Wong 
et  al.  53]  for  polygonal  (and  circular)  ducts  and  by  Wong  et  al.  [28] 
for  oval  ducts  with  different  long-short  axes  ratios.  The  critical 
and  neutral  radii  (rcr  and  re)  are  defined,  respectively,  as  the  outer 
radius  of  the  insulation  for  which  the  heat  transfer  rate  from  the 
circular  duct  is  a  maximum,  Qmax  (or  the  overall  thermal  resistance 
is  a  minimum,  £ftmin),  and  that  for  which  the  heat  transfer  rate 
from  the  circular  duct  is  the  same  as  that  of  a  non-insulated  duct 
(Qe  =  Qo)»  as  illustrated  in  Fig.  3.  Critical  heat  transfer  occurs  when 
the  outer  radius  of  an  insulated  circular  duct  is  less  than  its  critical 
insulation  radius.  This  phenomenon,  according  to  Wong  et  al.  [28], 
is  important  and  also  occurs  in  the  case  of  a  small-sized  insulated 
duct,  especially  in  conditions  of  low  ambient  air/gas  convective 
coefficients  (ft0).  The  study  showed  that  an  insulated  oval  duct  is 
similar  to  an  insulated  circular  duct  except  its  critical  and  neutral 
thicknesses  are  smaller  in  the  larger  oval  axes  ratio,  i.e.,  less 


Fig.  3.  The  critical  and  neutral  heat  transfer  phenomena  [48,3].  Qo,  heat  transfer 
rate  from  bare  pipe;  rcn  critical  radius;  re,  insulated  pipe  radius,  where  Q  is  the  same 
as  that  of  a  non-insulated  pipe;  rins,  insulated  pipe  radius;  and  £R,  overall  thermal 
resistance. 


insulation  material  would  be  required  if  the  long-short  axes  ratio 
of  the  oval  duct  increases  (i.e.,  for  a  flat-oval  duct).  Sahin  and 
Kalyon  [33]  specified  that  the  critical  radius  analysis  of  a  circular 
pipe  for  convection-type  heat  transfer  is  well  known,  and  the 
analysis  was  extended  to  radiation  heat  transfer  applications.  In 
that  study,  a  circular  pipe  in  which  the  outer  surface  is  subject  to 
both  convection  and  radiation  with  the  surroundings  was  handled, 
and  a  possible  critical  radius  for  various  parametric  conditions  was 
studied.  They  found  that  the  critical  radius  of  the  insulation  can  be 
calculated  for  a  circular  pipe  subjected  to  radiative  and  convective 
heat  transfer  environments  provided  that  certain  constraints  are 
satisfied  and  the  critical  radius  of  the  insulation  increases  when 
either  the  convection  or  radiation  decreases.  A  comprehensive 
study  on  a  critical  radius  analysis  of  pipes  and  ducts  was  carried 
out  by  Aziz  [35]  in  which  both  cylindrical  and  spherical  geometries 
were  studied.  In  addition,  cylindrical,  spherical,  rectangular,  equi¬ 
lateral  polygonal  and  eccentric  circular  insulation  geometries  were 
also  examined.  The  heat  transfer  characteristics  for  insulated 
and  non-insulated  spherical  containers  were  investigated  by  Wong 
et  al.  [56  .  By  both  considering  and  neglecting  the  effect  of  radiative 
heat  transfer  on  a  spherical  geometry  were  studied.  They  found  that 
the  effect  of  radiative  heat  transfer  cannot  be  ignored  in  conditions  of 
low  outside  convective  heat  transfer  coefficients  (such  as  ambient  air) 
and  high  surface  emissivities.  xins/r2  (insulation  thickness  divided  by 
outside  radius  of  sphere) =0.2  was  suggested  for  optimum  insulation 
thickness  in  practical  applications  because  of  that  the  difference  of 
heat  radiation  between  the  insulated  and  non-insulated  cases  was 
very  small  (within  3%)  for  between  xins/r2  =  0.2  and  xins/r2=l. 
Information  about  these  studies  and  the  particular  data  used  are 
summarized  in  Table  2. 

The  performance  of  an  insulation  material  is  mainly  deter¬ 
mined  by  its  thermal  conductivity,  which  is  dependent  on  the 
density,  porosity,  moisture  content,  and  mean  temperature  differ¬ 
ence  of  the  material.  Thermal  insulation  materials  in  various 
applications  are  exposed  to  significant  and  continuous  tempera¬ 
ture  variations  due  to  a  varying  ambient  air  temperature,  operat¬ 
ing  temperature  and  solar  radiation.  Therefore,  the  thermal 
conductivity  of  the  insulation  material  can  vary  due  to  changes 
in  both  moisture  content  and  temperature.  Condensation  in  the 
insulation  occurs  if  the  water  vapor  concentration  reaches  satura¬ 
tion.  Condensation  in  fibrous  insulation  is  a  serious  problem  in 
many  applications  such  as  building  insulation,  refrigerated  space 
envelopes  and  pipe/duct  insulation,  because  it  can  result  in  a 
drastic  reduction  in  the  thermal  insulation  of  the  system.  The  heat 


(/5 

0» 

•M 

O 

Z 


G 

o 

3 

C/5 

G 


C/5 

.G 

3 

n3 

S-h 

3 

u 


S-h 

u 

QJ 

sc 

-M 

WO 

c 

N 

>> 

3 

G 

03 

C/5 

_0J 

3 

G 

-i — * 
C/5 

QJ 

-G 


G 

O 

J3 

03 

G 

g 

03 

E 

3 

CN  G 

ii  £ 

3  3 

l-l 

H  CQ 


C/5 

a 

o 


•a 

a 

o 

u 

wo 

a 

'-3 

2 

01 

a 

o 


C/5 

a» 


at 


o 

at 

wo 

•a 

3 

a 

■w 

to 


C/5 

i-i 

01 

■M 

01 


2 

a- 


>5 

•a 

a 

■M 

C/5 


S-H 

<3  • 

^  C/5 

u  o 


^5  03 

’l-H  1-1 

Oj  s/5 


"O 

_g 

3 

G 


QJ 

S-H 

QJ 


C/5 

OJ 

S/5 

C/1 

CU 

G 

y 

U 


03 

S-H 

4-J 

G 

(U 

c 

"O 

G 

o3 

3 

u 


S-H 

u 

CU 
-G 
E— 1 


"O 

_g 

3 

G 


CU 

S-H 

(U 


s/5 

CU 

s/5 

C/5 

CU 

G 

y 

U 


03 

S-H 

■M 

G 

(U 

c 

T3 

G 

o3 

3 

u 


S-H 

u 

CU 
-G 
E— 1 


u 

03 

S-H 

03 

-G 

u 


CU 
X 
03 

-I— 1 

G 
(U 

<4-h 

lo 

G  T3 

2  Uh 
■M  O 

4- . 

03  -Q 
^  OJ 

»T“  +_l 

^  03 

3  .£P 

(U  ■ 

E 

5- H 

,o 


s/5 
CU 
> 
G 

M— I  •  >-H 

J-H 

CU  % 

a  2 

■4—4  > 

o  ^ 
c  tl 

£  ^ 
G  5 
.2  o 

-M 

03  03 
N  <u 

I  3 

■D  G 
cx  £ 

o  .2 


cu 

S-H 

(U 


ttJ 

n. 

J- 


u 

G 

03 

G 

O 

u 

cu 

a 


03 

c 

03 


G 

2 

u 

C-j-H 

CU 

o 

u 

G 

_g 

'■D 

u 

cu 

C  3 

o  oj 
u  u 

S-H  CU 
CU  ' — 1 

C  W> 

p  CU 

3  c 


03 

G 

03 

S'  ^ 

Vs  X 

| 

m  ^  > 

ro  H 

°?  -G° 

in  ^ 

II  ^ 

-3  E 

£ 

G  o 

(U 

u  I 
UG  ^ 


s/5 


CU 

O 

U 

G 

_g 

’•M 

u 

(U 

> 

G 

o 

u 

S-H 

(U 

G 

C 


cu 

.SO 

J1 


s/5 

CU 


U 

G 

03 

G 

O 

u 


ro 

od 


o 

-e 

(N 


LD 

O 


-c 

u 

o 

o 


h? 

u 

o 

o 

o 


>  cu 

>  CX 
o  -s 
I-J  cx 


cu 

S-H 

03 

3 

cr 

S/5 

s-T 

J3 

3 

bJO 

C 

.2  c/ 

£  s. 

•s  'a 

G  'T 

O  u 
bO  ■M 

i?  ^ 

g  3 
ft  a 
o 

03  bO 
Th  03 
G  -M 
bJO  G 
<U  OJ 

(xs  a 


cu 

S-H 

03 

G 

cr 

C/5 

C/1 

s-T  CU 

-2  9? 
g  'S, 

WO  Lh 

5  * 
■c  => 
E 

3 

G  03 
O  G 
bO  03 
— N 

g  3 
cx  a 
o 

03  bJQ 
—  03 

G  -M 
WO  G 
CU  cu 
Dei  a 


>s 

S-H 

03 

> 

c/i 

_g 

'•D 

03 

S-H 

C/5 

CU 

X 

03 

-M 

S-H 

O 

J3 

‘f  CD 

WO  03 

2  c 
O  o3 


U 

G 

03 


G 

CU 

(U 


G  -M 
>  CU 

O  J3 


CD 

05 

ro 

I 

ro 

T - 1 

CN 


3 


of 

O 

o' 


r^ 


m 

ro 

O 

o' 


to 

E 


to 


o 

3  jH 

■a  -5 

-is? 


3  ^ 

O  S 
u 


WO 

G  , — , 
o  00 

^  21 


WO  03  OT  WO  G 

r-  c^s  i— 


CU 


i£L  § 


cu 


0.  Kaynakli  /  Renewable  and  Sustainable  Energy  Reviews  30  (2014)  184-194 


191 


03 

03  <y 

CL) 


D 

Uh 

c2 

(U 

> 

'■M 

u 

cu 

> 


-2  o 


S-H 

CU 


c/i 

_G 

3 

03 

S-H 


G  _ 

O  a3 
U  u 

o  '£ 

4- 1  S-H 

c  u 

O  03 
•£3  G 

3  ^ 

03  03 
03  (U 

S-H 

G  (U 
•  — i  03  -  • 

5- 1  3 

G 
O 
u 


.CU 


c/i 

G 

03 

S-H 


3  QJ 

3 

G 


03 

QJ 

-G 

QJ 

> 

'-M 

.2 

3 

03 

CP 


c r\ 

C/5  C/5 

S  I 

CX  c- 

S-H  1—1 

03  P 


G 

u 

S-H 


C/1 

G 

03 

S-H 


CN 

o' 


cu 

to 


I 

CN 

o' 


o 

-C 


ro 

o' 

l 

O 


5 


C/l 

-l-l 

U 

G 

03 

S-H 

-2 

3 

a 

S-H 

u 


S-H 

p 

C-t-H 

C/l 

G 

03 

S-H 


03 

QJ 

-G 

QJ 

.> 

'-D 

.2 

3 

03 

S-H 

03 

G 

03 

QJ 

U 

S-H 

3 

03 

G 

o3 


2  T3 

G  E 

3  oj 
c  3 

^ ^  C/l 

C  C 

o  o 
3  u 


QJ 

S-H 

QJ 


U 
QJ 

g 

s  ^3 
3  ^ 

o 

CO 


o 

E 


S-H 

,o 


S-H 

,o 


CN 

CO 

05 


•E  E 

rs  cd 

E—<  of 

t  00 
S  II 

O  o 
jc 

3  -a 
II  c 

o  n3 

*a  (n 

c  E 


03  _ 

N  ^ 


(N 


CN 

CN 


CN 

uo 

oo 

[ 

o 

rd 


o 

u 

o 

o 

00 


eE1 


o 

*C 

c” 

_g 

*-M 

u 

QJ 

> 

G 

O 

u 

3 

S-H 

G 

-i — * 

03 


-  S-H 

3  3 

C  G 

o  y 

WO  .b 
>>  u 

3  cj 

CX  -C 


03 

S-H 

QJ 


G 

QJ 
U 

3  QJ 

’3  3 
a*  a 

QJ  03 

_ -  S-H 

3  3 
3  a 

s-h  bJJ 
03  G 
G  03 

•  —  -i— i 

u 
QJ 

S-H 


03 

QJ 

S-H 

QJ 

3 

'c/5 

G 

O 

u 

C/5 

03 


S-H 

P 

C-t-H 

C/l 

G 

03 

S-H 


03 

QJ 

-G 

QJ 

> 

'-M 

.2 

3 

03 

S-H 

^>5 

3 

o 


o 

-G 

S-H 

.O 


C/l 

G 

O 


P 

G 

O 

u 


^  o 


O 
U 

S-H 

,P 

o 

-  CN 

E 


-e  ^ 
u 


CN  > 

ro  > 
I  ro 

ro  °°' 


CT5 

o' 

I 

O 

II 

w 

c” 

o 


t> 

3  ^ 
roo  - 

u 

„  o 

u  o 

°  CN 
O  | 
ro  1 


o  u 

u  o 
P  o 

QJ 

u 

S-H 


c ri 

G 

O 


P 

G 


O  /-  O 

i  h-i  D 


U 


c/i 

OJ 

’C 

C/l  -4— » 

cu  cu 
^  p 
S  o 

-  & 

G  _H 
O  03 

■43  3 

03  s-h 

G  ^ 

G  ft 
.G  on 


E 

3-  CD 


3  II 

O'  E 

II 

cn  » 

>>  H 
3rt 

rb  -s-t 
S-h  C/5 
OJ  >, 

G 


E 

u 

of 

I 

m 

CN 

o' 


IN  c 

II  E 

tu  _ _ _ 

.a,  O 
n.  C75 


^  3 


s-?  m 


QJ 


E 

E 


QJ 

U 


ro  c/ 

2°.  c  o 
g  o  o  o 
3  II  p  cn 

P  1 '  S— t  II 

G  E  *03 

£ 


W)  G  ^ 

^  LO 

p  D  Hi 


and  moisture  transfer  characteristics  in  fibrous  insulation  were 
investigated  by  Leskovsek  and  Medved  [57  ,  Fan  and  Wen  58]  and 
Choudhary  et  al.  [59].  Leskovsek  and  Medved  [57  found  that  a 
relatively  small  mass  of  water  in  the  insulation  matrix  can  result  in 
a  significantly  increased  average  heat  flux.  Moreover,  Fan  and  Wen 
[58]  found  that  the  initial  water  content,  the  thickness  of  the 
fibrous  insulation  and  the  environmental  temperature  are  the 
three  most  important  factors  influencing  the  heat  flux.  Choudhary 
et  al.  [59]  studied  the  condensation  around  a  cold  pipe  and 
the  removal  of  the  condensate  by  a  hydrophilic  wick  fabric  with 
downward  orientation  by  the  combined  effect  of  capillary  and 
gravity  actions.  As  a  result,  the  presence  of  the  wick  appears  to 
significantly  reduce  the  amount  of  liquid  water  in  the  insulation 
and  to  maintain  the  conductivity  of  the  insulation. 


4.  Economic  analysis 


In  the  literature,  several  economic  analysis  methods  are  used  to 
optimize  the  thermal  insulation  thickness  for  different  applica¬ 
tions,  such  as  walls,  pipes  and  ducts,  and  to  calculate  the  payback 
period.  One  of  the  economic  methods  is  the  Simple  Payback 
Analysis.  This  method  is  based  on  the  time  required  to  repay  the 
initial  capital  investment  with  the  operating  savings  attributed  to 
that  investment.  The  main  drawback  of  this  simple  analysis  is  that 
it  does  not  take  into  account  the  inflation  and  the  time  value  of 
money,  which  are  very  important  financial  considerations  [60,61]. 
An  economic  assessment  is  generally  performed  using  a  life-cycle 
cost  (LCC)  method  during  each  period  of  time.  In  this  analysis,  the 
amount  of  energy  savings  due  to  the  insulation  over  a  lifetime  is 
evaluated  in  the  present  value  using  the  present  worth  factor 
(PWF).  The  PWF ,  which  depends  on  the  inflation  rate  (z)  and  the 
discount  rate  (d),  is  calculated  as  follows  [18,25,62]: 


PWF  = 


(if  d  i-  i ) 


(23) 


PWF  = 


LT 
1  +d 


(if  d  =  z) 


(24) 


The  LT  could  be  assumed  to  be  10  years  [25,26,37,38],  15  years 
[24]  or  20  years  [36]. 

The  P\-P2  method  presented  in  Duffie  and  Beckman  [63]  is 
used  in  several  studies  [26,37,38  .  This  method  simplifies  the 
economic  analysis  by  concentrating  all  of  the  economic  para¬ 
meters  into  two  parameters,  and  P2.  Pi  and  P2  are  determined 
from  the  economic  parameters  (e.g.,  interest  rate,  inflation  rate, 
period  of  economic  analysis).  Pt  is  the  life  cycle  energy  related  to 
the  market  discount  rate  d  (for  the  value  of  money),  the  inflation 
rate  z  (for  the  energy  cost),  and  the  economic  analysis  period  (or 
the  technical  lifetime  of  the  applied  insulation  in  years  LT  [26,64]). 
The  value  of  Pi  can  be  calculated  by  the  following  equation: 


lt  (i+jy-i 
j= i  (i+dy 


i 

(d  -  i ) 

LT 

1+2 


z  ^  d 


z  =  d 


(25) 


Typically,  for  a  period  of  economic  analysis  of  LT  years,  the 
value  of  Pi  varies  between  LT/2  and  LT.  P2  is  the  ratio  of  the  life 
cycle  expenditures  incurred  due  to  the  additional  capital  invest¬ 
ment  to  the  initial  investment,  which  can  be  calculated  by  the 
following  equation: 

P2  =  DP  +  (  1  —  DP)P,  +MSP,  -  Rv  (26) 

(1  +d) 

where  DP  is  the  ratio  of  the  down  payment  to  the  initial  investment, 
Ms  is  the  ratio  of  first  year  miscellaneous  costs  (maintenance, 
insurance,  and  other  incidental  costs)  to  the  initial  investment,  and 
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Rv  is  the  ratio  of  the  resale  value  at  the  end  of  the  economic  period  to 
the  initial  investment.  Typically,  the  value  of  P2  varies  between 
0.5  and  1.0.  If  no  additional  capital  is  invested  other  than  the  initial 
investment  (such  as  maintenance  and  operation  costs),  P2  can  be  set 

to  1.0  [17,37,38,64]. 

It  is  noted  that  because  the  chosen  economic  parameters  can 
strongly  affect  the  present  worth  of  the  total  life  cycle  cost  (and 
thus,  the  optimum  thermal  insulation  thickness),  they  must  be 
accurately  determined  [47]. 


5.  Investigation  of  parameters  affecting  optimum  thickness 

A  practical  application  example  for  optimizing  the  thermal 
insulation  thickness  on  a  pipe  is  conducted  using  the  Pi-P2 
method,  and  the  factors  affecting  the  economic  thickness  are 
investigated  in  the  present  study.  In  the  present  analysis,  a  circular 
pipe  subjected  to  convective  heat  transfer  to  the  surrounding 
environment  is  studied.  It  is  assumed  that  the  carrier  pipe  has  a 
100  mm  nominal  diameter  (112  mm  outer  diameter),  is  insulated 
with  mineral  wool  with  a  thermal  conductivity  of  0.040  W/mK 
and  operates  at  120  °C  for  3000  h  per  year.  The  useful  lifetime 
of  the  system  is  assumed  to  be  10  years  [25,26,37,38],  and  the 
inflation  and  discount  rates  are  4%  and  5%,  respectively.  The 
convective  heat  transfer  coefficient  at  the  surface  of  the  insulation 
is  set  to  a  typical  value  of  10  W/m2I<  [29,34,47  .  The  particular 
parameters  used  in  the  calculations  are  presented  in  Table  3.  As 
the  insulation  thickness  on  a  pipe  increases,  the  heat  loss  per  unit 
length  of  pipe  decreases  (and  thus,  the  energy  requirements),  as 
expected.  The  variation  in  the  cost  curves  (energy,  insulation  and 
total)  with  the  thermal  insulation  thickness  and  the  optimum 
economic  thickness  of  the  insulation  can  be  observed  in  Fig.  4. 
With  an  increasing  insulation  thickness,  the  total  life  cycle  energy 
cost  decreases.  In  contrast  to  flat  surfaces,  such  as  the  external 
walls  of  a  building,  the  increase  is  not  linear  because  of  the 
amount  (volume)  of  insulation  used  on  the  pipe.  The  optimum 
insulation  thickness  that  minimizes  the  total  life  cycle  cost  over  a 
lifetime  of  10  years  is  found  to  be  0.1136  m. 

The  variations  in  the  optimum  insulation  thickness  with  the 
convective  heat  transfer  coefficients  inside  and  outside  of  the  pipe 
are  shown  in  Figs.  5  and  6,  respectively.  Typical  value  ranges  of  the 


Table  3 

The  parameters  used  in  the  calculations. 


Parameter 

Value 

Pipe 

Inner  radius,  rA 

0.05  m 

Wall  thickness,  t 

0.006  cm 

Outer  radius,  r2 

0.056  m 

Thermal  conductivity,  kpipe 

54  W/mK  [36] 

Inner  heat  transfer  coefficient,  hi 

100  W/m2I<  [34] 

Outer  heat  transfer  coefficient,  ha 

10  W/m2I<  [29,34,47] 

Working  fluid  temperature  in  pipe,  Tf 

120  °C 

Outside  air  temperature,  T0 

15  °C 

Energy  source  (natural  gas) 

Lower  heating  value,  Hu 

34526  kj/m3 

Cost,  Cen 

0.4  US$/m3 

System  efficiency,  rj 

0.93 

Thermal  insulation  (mineral  wool) 

Thermal  conductivity,  kins 

0.04  W/mK  [37,47] 

Cost,  Cins 

140  US$/m3 

Economic  parameters 

Inflation  rate,  i 

0.04 

Discount  rate,  d 

0.05 

Lifetime,  LT 

10  [25,26,37,38] 

Operating  time,  N 

3000  h 

Fig.  4.  The  variation  in  life  cycle  costs  with  insulation  thickness. 


Fig.  5.  The  optimum  insulation  thickness  vs.  the  inside  heat  transfer  coefficient. 

forced  convective  heat  transfer  coefficient  inside  the  pipe  for  air 
or  water  flow  are  also  given  in  Fig.  5.  Furthermore,  for  a  pipe 
surrounded  by  air,  typical  convective  heat  transfer  coefficients  at 
the  outer  surface  of  the  insulation  for  free  and  forced  modes, 
which  are  related  to  the  wind  velocity,  are  also  shown  on  Fig.  6.  As 
shown  in  the  figures,  the  optimum  insulation  thickness  increases 
with  increasing  heat  transfer  coefficients.  Insulation  provides  the 
primary  thermal  resistance  against  heat  loss.  Therefore,  if  the  heat 
transfer  coefficient  increases,  the  total  thermal  resistance  of  the 
pipe  decreases,  and  this  results  in  an  increased  insulation  thick¬ 
ness.  The  variation  in  the  insulation  thickness  is  higher  at  low  heat 
transfer  coefficients.  Because  the  internal  heat  transfer  coefficient 
is  influenced  by  the  fluid  properties,  the  type  of  working  fluid  in 
the  pipe  (air,  water,  etc.)  is  essential.  Moreover,  because  low  heat 
transfer  coefficients  inside  the  pipe  or  the  duct  are  generally 
encountered  for  gas  flows  such  as  air,  the  flow  velocity  of  the 
gas  fluid  in  a  pipe  or  a  duct  should  be  determined.  In  general,  the 
effect  of  the  inner  heat  transfer  coefficient  on  the  optimum 
insulation  thickness  is  slightly  lower  than  that  of  the  outer  heat 
transfer  coefficient,  as  seen  in  Figs.  5  and  6. 

The  working  fluid  temperature  in  a  pipe  is  one  of  the  important 
factors  affecting  the  required  insulation  thickness.  Fig.  7  presents 
the  effect  of  the  fluid  temperature  on  the  optimum  insulation 
thickness  for  different  thermal  conductivities  of  insulation  materi¬ 
als.  In  the  calculations,  the  values  of  the  other  parameters  were 
taken  as  given  in  Table  3.  Because  an  increase  in  the  fluid 
temperature  and  in  the  conductivity  of  the  insulation  leads  to  an 
increase  in  the  heat  loss  rate  from  the  pipe,  a  greater  amount  of 
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Fig.  6.  The  optimum  insulation  thickness  vs.  the  outside  heat  transfer  coefficient. 


Fig.  7.  The  optimum  insulation  thickness  vs.  the  fluid  temperature  in  a  pipe. 


Fig.  8.  The  optimum  insulation  thickness  vs.  lifetime. 


insulation  is  required.  Thus,  the  economic  insulation  thickness 
increases  with  an  increase  in  the  fluid  temperature  and  in  the 
thermal  conductivity  of  insulation. 

The  expected  useful  lifetime  of  a  system  and  the  operating  time 
in  a  year  are  among  the  important  factors  in  an  economic  analysis. 
In  the  literature,  the  lifetime  was  generally  assumed  to  be  10 
years;  however,  various  periods  ranging  from  5  years  to  30  years 
were  also  used  [19,20,24,25,36].  To  investigate  the  influence  of  the 
operating  time  in  a  year,  a  wide  range  from  1000  h  to  8760  h, 


Fig.  9.  The  optimum  insulation  thickness  vs.  pipe  diameter  (D2). 

which  covers  the  entire  year,  is  considered.  Fig.  8  shows  that  the 
optimum  insulation  thickness  increases  with  an  increase  in  both 
the  lifetime  and  the  annual  operating  hours.  The  increases  are  not 
linear;  therefore,  as  the  lifetime  or  the  operating  time  increases 
by  three  times  (from  10  years  to  30  years  or  from  1000  h  to 
3000  h),  the  optimum  thickness  does  not  increase  in  the  same 
ratio.  Moreover,  the  annual  operating  time  is  slightly  more 
effective  than  the  lifetime  on  the  optimum  insulation  thickness. 
The  effect  of  pipe's  outer  diameter  on  the  optimum  insulation 
thickness  is  shown  in  Fig.  9.  In  the  calculations,  the  values  of  the 
parameters  were  taken  as  given  in  Table  3,  except  for  the  value  of 
the  pipe  diameter.  Because  the  heat  transfer  surface  area  increases 
with  an  increasing  pipe  diameter,  the  required  insulation  thick¬ 
ness  increases  to  limit  the  heat  transfer  rate.  The  increment  has  a 
decreasing  slope,  as  shown  in  Fig.  9.  Therefore,  as  the  pipe 
diameter  increases,  a  relatively  thinner  insulation  thickness  is 
required. 

6.  Conclusions 

In  recent  years,  a  significant  increase  in  the  number  of  studies 
that  have  been  done  on  the  performance  of  thermal  insulation 
systems  applied  to  flat  surfaces,  cylindrical  surfaces  and  different 
geometries  has  been  observed.  Reducing  the  energy  consumption 
in  industry  is  important  when  considering  the  limited  energy 
resources  and  the  environmental  concerns.  This  study  presented  a 
literature  review  on  the  thermal  insulation  applications  in  pipes 
and  ducts,  and  focused  particularly  on  the  determination  of  an 
optimum  insulation  thickness  considering  economic  criteria.  The 
insulation  materials,  the  fuel  types,  the  pipe/duct  geometries,  the 
fluids  types  and  the  operating  conditions  encountered  in  the 
literature  were  compiled  and  presented.  The  effects  of  thermal 
insulation  and  other  design  and  operating  parameters  on  the 
energy  loss  from  pipes  and  ducts  were  investigated.  This  study 
could  serve  as  a  useful  source  of  information  for  researchers 
because  it  includes  a  review  of  the  thermal  insulation  applications 
particularly  in  pipe  and  duct  systems,  the  optimization  procedure 
for  determining  the  economic  insulation  thickness  and  a  practical 
application  example. 
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